Phase diagram and influence of defects in the double perovskites 
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The phase diagram of the double perovskites of the type S^-^LaxFeMoOe is analyzed, with and 
without disorder due to antisites. In addition to an homogeneous half metallic ferrimagnetic phase 
in the absence of doping and disorder, we find antiferromagnetic phases at large dopings, and other 
ferrimagnetic phases with lower saturation magnetization, in the presence of disorder. 
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Introduction. The, double perovskite Sr 2 FeMo06 
and related materials!!] are good candidates for mag- 
netic devices, as they combine a high Curie temper- 
aturemand a fully polarized (half metallic) conduction 
bandfl A^s^eaL-tJieise .materials are being extensively 
studiedMMffl'§BNolB.Byy 

The magnetism of these compounds arises from the 
Fe 3+ , S — 5/2 core spin, while the charge state of the 
Mo ion is 5+. Spatially, the Mo and Fe ions occupy 
two interleaving FCC lattices (sodium chloride struc- 
ture). The conduction band contains one electron per 
unit cell, which tends to be antiparallel to the Fe spin. 
Experiments suggest that, in many samples, the satu- 
ration magnetization is less than the expected 4//g per 
formula unit. This effecLisjjsu^l^^pdb^d to the pres- 
ence of antisite defects^M-BBEMiSU where, due to 
the similarity of their atomic radii, Mo ions are randomly 
placed on the Fe sublattice and conversely. Notice that 
when a Fe ion is misplaced, with high probability it will 
have a Fe ion among its first neighbors, enhancing direct 
antiferromagnetic (AFM) superexchange with respect to 
the ideal structure. The strength of this coupling can be 
inferred from the compound LaFeC-3, which has the same 
structure, but where the Mo ions have been substituted 
by Fe 3+ . LaFe0 3 is known to be AFMjl3 with a Neel 
temperature of Tn = 720K. 

The Sr ions in Sr 2 FeMoOg can be substituted for triva- 



lent cations, like La, leading to Sr 2 _ a; La a ;FeMo06J 
These compounds have 1 + x electrons per formula unit 
in the conduction band. These doped materials tend 
to have a higher Curie temperature. Notice that one 
can also consider the substitution with a monovalent ion 
(i.e. Sr— > K, S^-^K^FeMoOe ), which takes one elec- 
tron from the conduction band, leaving 1 — x electrons 
per formula unit. Hence in this paper negative x will 
actually refer to substitution with a monovalent ion. 

The model. Band structure calculations have shown 
that the conduction band can be described ir^-terms of 
hybridized f 2g orbitals at the Mo and Fe sitesBo If one 
considers the i 2ff orbitals of both spin orientations at the 
Fe sites, the model leads to a highly correlated system, 
where an on site Hund's coupling and a Hubbard repul- 
sive term have to be addedE^uil'L 2 ] In the following, we 
will consider the magnetic phase diagram only, and ne- 
glect the possible existence of a metal- insulator transition 
when the ratio between thej,bandwidth and the Coulomb 
term is sufficiently small.EilEa We consider that the con- 
duction band is built up of the three i 2s orbitals at the Fe 
sites with spins oriented antiparallely to the Fe moment, 
and the six i 2g orbitals at the Mo sites (see below). 

We denote the destruction operator on xy orbitals with 
spin + or — at lattice site r as F xy -±. r , M xy -±. r (F for 
Fe and M for Mo), and so on. The spin and number 
operators on a given Fe site are: 



a,/3=± 

N r = ^""^ {Fly;a; r FxVWr + ^xz-.a :r ^xz-.a; r + Fy Z . a . r F yZ]a]r ) . 
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a=± 



Analogous definitions hold for the Mo atoms. Given the large spin value (S=5/2) of the the localized Fe core spins, 
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we treat them as classical, with polar coordinates 

4> = (sin 9 cos ip, sin sin ip, cos 8) (3) 

As mentioned above, we only consider the Fe orbitals 
with spin antiparallel to <f>, which amounts to assume 
that the Hund's coupling at the Fe ions is much larger 
than the other interactions. Thus, we define up and down 
orbitals, /i and fa, with respect to the local 5/2 spin: 

6 
F + = cos -fx + sin -f 2 , 

F- = sin l^ v fi- cos ^>/ 2 , (4) 

and we neglect all terms including the /i operators. For 
the sake of brevity, in the following we set f% = /. Then, 
the Hamiltonian, in the absence of disorder and neglect- 
ing direct hopping terms between Mo orbitals (see be- 
low), can be written as: 

H = IC xy +IC yz +IC xz -fi J2 M* e -(»+A)J2 (5) 

■ncvvn 7* odd 

with 

Q 

fcxy = ^Mo-Fc (sin fl y . r M xy . + . r+ u + h.c.) ... 

T>£Fe lattice " 
u = e x ,e y 

(6) 

where, for brevity, we omit the analogous hoppings from 
the sites belonging to the Mo sublattice to the Fe sites. 
Analogous expressions are found for the kinetic energy on 
the xz and yz planes. Finally, we add a direct hopping 
between Mo orbitals. These hopping terms fldwhgjfe to 
three separate two dimensional Hamiltonians.E3'EilE3 The 
substitution of Mo ions for Fe ions leads to direct Fe-Fe 
hopping terms, ( we take ^fc-Fo = ^Mo-Fc), and also to 
the inclusion of an AFM exchange, Jfc-Fc- Thus, the 
model is defined by the parameters iMo-Fc, tivio-Mo, A, /i 
and Jfc-Fc- There are nine orbitals per unit cell, three 
at the Fe sites, and six at the Mo sites. 

The occupancy of the conduction band depends on 
the value of the chemical potential, fi, and it varies 
from one electron to two electrons per unit cell in 
Sr2- x La x FeMo06, < x < 1. We neglect interactions of 
the electrons within this band (as discussed below, the 
number of electrons at the Fe sites is always less than 
one). 

We will use iFc-Mo as our unit of energy (iFe-Mo ~ 
0.35eV from band structure calculations). We take 
^Mo-MoAfc-Mo = 0.25, A = and Jfc-FcAfc-Mo = 0.1. 
The value of A implies a relatively large hybridization 
of the Fe and Mo orbitals— which seems consistent with 
Hartree-Fock calculations.Ej Jf c -Fc is chosen so as to re- 
produce the Neel temperature of LaFeOs. We have not 
made a comprehensive study of the dependence of the 
results on the tight binding parameters, but the calcula- 
tions made so far indicate that the qualitative features of 



the phase diagrams to be discussed below are not strongly 
dependent on the choice of parameters. In the absence 
of disorder, this model is basically equivalent to the one 
studied by Chattopadhyay and MillisEj in the context 
of Dynamical Mean Field Theory, although we shall use 
Variational Mean Field (see [2^ for a comparison between 
the two methods). T\]&Mi^hL^^^sJm our considering 
of the disorder effectslfflflll3'BBEJO: with probability 
y we misplace an Fe ion onto the Mo sublattice (and 
conversely) without any spatial correlations (y is just the 
antisite density). It is clear that y = 0.5 corresponds to 
full disorder on the location of the Fe and Mo ions, while 
y > 0.5 is equivalent to 1 — y with the Fe and Mo sublat- 
tices interchanged. Vacancies can be equally considered, 
but explicit calculations showed that they have a much 
milder effect on the phase diagram. 

Method of calculation. We use the method developed 
for double exchange systems in Ref. We assume that 
the Fe core spins are classical. At a given temperature, 
we average over spin configurations obtained by assum- 
ing that there is a magnetic field acting on the spins. 
The magnitude of these fields are variational parameters, 
which are taken so as to minimize the free energy. Given 
a spin configuration, the electronic states are calculated 
exactly, and the electronic contribution to the free energy 
is obtained by integrating the density of states. As the 
Fe spins are distributed in a three dimensional lattice, 
and the electrons lead to effective interactions with the 
cubic symmetry, we think that our mean field ansatz for 
the spin configurations is sufficient. This method is in 
excellent agreement with more precise Mente Carlo cal- 
culations for the double exchange modeled We solve the 
Hamiltonian in lattices with up to 512 x 512 x 512 sites 
(note that the calculation of the electronic wave functions 
requires only the diagonalization of the Hamiltonian in 
a 512 x 512 square). For these sizes, the disorder due to 
antisites is self-averaging. 

The adequacy of our technique depends on the ansatze 
made for the possible spin configurations. We have con- 
sidered four possible phases: i) the paramagnetic (PM) 
phase, ii) the ferrimagnetic (FI) phase, where all Fe spins 
are parallel, and the spins of the electrons in the conduc- 
tion band are antiparallel to the Fe spins, iii) an AFM 
phase, where the Fe spins in neighboring (1,1,1) planes 
are antiparallel, and iv) a different ferrimagnetic (FIP) 
phase where the Fe spins are aligned ferromagnetically if 
the Fe are in the correct positions, and antiferromagneti- 
cally if the Fe ions occupy Mo sites because of the antisite 
defects. In the absence of disorder, we have checked that 
other phases with canted spins have higher free energy. 
Note that the above ansatze define the average magneti- 
zation at the Fe sites, but that thermal fluctuations are 
also included. 

Results. The phase diagram of Sr 2 _ a; La 2 ;FeMo06, as 
function of x and temperature, is shown in Fig. [I] for 
different concentrations of antisites. 

In the absence of defects, we find that T c decreases 
with increasing doping of the conduction band, in agree- 
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antisites: 0% - 




Figure 1: Phase diagram of S^-^LaxFeMoOe as function of 
x and temperature for different concentrations of antisite de- 
fects. Negative x actually means Sr 2 _| I |K| :E |FeMo06. In both 
cases, the density of carriers in the conduction band is 1 + x. 
Phase-separation regions are found between the FI and AFM 
phases (upper panel), and between the FIP and AFM phases 
(middle and lower panel). 



ment with Ref. 2^. At high, but still reasonable, dopings 
we find the ordered AFM phase described above. The 
phase transitions are first order, with regions of phase 
separation between them. For x ~ 0, the spins of the 
electrons at the Mo orbitals are antiparallel to the Fe 
core spins. We ascribe the tendency toward phases with 
zero magnetization, upon increasing doping, to the oc- 
cupancy of the Mo orbitals which are aligned parallel to 
the Fe spins. 

The presence of antisite defects changes significantly 
the phase diagram: i) The FI phase is replaced by the 
FIP phase, where the spins at the Fe sites at the de- 
fects are antiparallel to the overall magnetization, ii) the 
ordered AFM phase is strongly suppressed, and iii) the 
value of Xa increases as the concentration of antisites also 
increasesEa, iv) the dependence of 7c with the number of 
electrons in the conduction band is more pronounced in 
the presence of antisites. 

These effects are associated to the direct AFM interac- 
tion between spins at Fe ions which are nearest neighbors. 
These interactions play no role in perfect materials. The 
antiferromagntic interaction can be easily shown to be 
equivalent to a ferromagnetic one for the atoms in the Fe 
sublattice. Thus, superexchange enhances the tendency 
toward a ferromagnetic order in the original Fe sublat- 
tice. This effect is independent of the number of electrons 
in the conduction band. The saturation magnetization, 
on the other hand, is reduced. 

Fig. H gives the occupancies of the different orbitals as 
the number of electrons in the conduction band is varied. 
Most of the charge is in the Mo orbitals. The variation 
is not-linear, indicating that a rigid band picture is not 
valid.O There are sharp changes at the phase transitions. 



Figure 2: Occupation of the Mo-f, Moj, and Fej, as function 
of the doping of the conduction band. The curves give the 
occupancies for a 10% density of antisites defects. Note that 
in phases with no net magnetization, the occupancies of the 
Mo T and Mo; levels are the same. 
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Figure 3: Low temperature magnetization (/ib per formula 
unit) of Sr2FeMoC>6 as function of the concentration of anti- 
site defects, y. Experimental results are from Refs. B H , E3 . 



At low temperatures, the spins at antisites tend to be 
antiparallel to the magnetization, as shown in Fig. []]. 
This implies that the saturation magnetization is reduced 
with respect to the ordered case. The total magnetiza- 
tion of the core spins and the conduction electrons, is 
shown in Fig. |[ The calculated magnetization is well 
fitted by the line Ms — (4.0 — 7.7y)jUB, where y is the 
antisite density. Experimental results from Refs. P, Oll3 . 
are added for comparison. Note that the decrease in the 
magnetization does not lead to a lowering of the Curie 
temperature, as discussed above. 

Conclusions. We have studied the magnetic phase dia- 
gram of the doped double perovskites, S^-zLa^FeMoOe. 
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We have analyzed the influence of antisite defects, on the 
phase diagram. 

In clean systems, we find that, as the number of 
electrons in the conduction band increases, the criti- 
cal temperature decreases, in agreement with previous 
calculations EB This variation is due to the increased fill- 
ing of the Mo-f band, which reduces the double exchange- 
like mechanism which tends to align the Fe moments. At 
sufficiently high dopings, we find ordered phases with- 
out net magnetization, which enhance the derealization 
of both the Mo-|- and Moj bands. The transitions be- 
tween these phases tend to be first order, with regions 
of phase separation between them. Electrostatic effects 
will prevent the existence of phase separation at macro- 
scopic-spales, leading to a domain structure at mesoscopic 
scalescS. 

Antisite disorder induces significant changes in the 
phase diagram. The ordered ferrimagnetic phase is re- 
placed by a different ferrimagnetic phase where the Fe 
spins at defects are antiparallel to the bulk magnetiza- 
tion (the FIP phase, see Fig.|f|). Antiferromagnetism at 
finite dopings is suppressed. The saturation magneti- 
zation in the FIP phase is reduced, although the Curie 
temperature tends to increase with the number of Fe in 
Mo positions, due to the direct AFM exchange between 
Fe ions which are nearest neighborsEa. 

Note that, in order to study compounds with different 
number of carriers in the conduction band, the presence 
of vacancies and changes in the Fe - O - Mo bond angles, 
Fe/Mo - O distance, and in the energy splitting A can 



influence the results. These effects need to be extracted 
from the available experimental data and incorporated in 
the model Hamiltonian, eq. (||). 

We have not studied transport properties, although 
it seems likely that the variation of the magnetic struc- 
ture near defects-will lead to significant changes in a half 
metallic system.113 We have also not analyzed other ef- 
fects of the electron-electron interaction, such as the exis- 
tence of a Mott transition to an insulating state, found in 
the related compound Sr2FeW06-c2l We think, however, 
that our model includes all relevant interactions required 
to study the magnetic properties of the metallic state of 
double perovskites. Similar models provide a good under- 
standing of the magnetic properties of the half-metallic 
manganite oxides (such as Lai-^Ca^MnOa)]^'^ 

In summary, we find a rich phase diagram for 
Sr2- x Laa;FeMo06, which is significantly modified in the 
presence of defects. Our results seem consistent with ex- 
isting experimental data. 
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